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ABSTRACT 


The water deer, Hydropotes inermis (Cervidae), is native to China and Korea and has two subspecies of the Chi- 
nese water deer (Hydropotes inermis inermis) and Korean water deer (Hydropotes inermis argyropus). To date, 
only the Korean water deer has been reported in South Korea. In this study, however, an intraspecific phylogeny 
and haplotype analysis based on mitochondrial cytochrome oxidase I indicated that both Korean and Chinese 
water deer are found in South Korea. The populations of the two Korean genetic lineages did not show distinct 
geographic distributions. Further morphological studies on the Korean water deer will be required to confirm its 


taxonomic status. 
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INTRODUCTION 


The water deer, Hydropotes inermis (Cervidae), is native to 
China and Korea and has two of subspecies of the Chinese 
water deer (Hydropotes inermis inermis) and Korean water 
deer (Hydropotes inermis argyropus) (Grubb, 2005; Harris 
and Duckworth, 2008). The Chinese water deer is distribut- 
ed to only the lower Yangtze Basin (Xu et al., 1996; Wang, 
1998; Hu et al., 2006), and it is listed as a vulnerable species 
in the IUCN Red List of Threatened Species (Harris and 
Duckworth, 2008) and the China Red Data Book of Endan- 
gered Animals (Wang, 1998). In contrast, the Korean water 
deer is widely distributed throughout the forests of South 
Korea, and it has been designated as a wildlife pest by the 
Korean Ministry of Environment. 

According to morphological studies (Allen, 1940; Won, 
1967; Won and Smith, 1999), only the Korean subspecies 
has been recognized in South Korea to date. However, a 
recent analysis of the mitochondrial Cyt b gene and its con- 
trol region (Koh et al., 2009) suggested the presence of two 
genetic lineages in Korean water deer populations, and one 
of these lineages was grouped with the Chinese subspecies. 

The objective of the present study was to confirm the pre- 
sence of the two genetic lineages in the Korean populations 


by analyzing mitochondrial cytochrome oxidase I (COI) 
gene sequences. We analyzed the intraspecific relationships 
in Korean water deer by using tissues collected from 18 
road-killed individuals from two northern forests (SC: Sok- 
cho, YG: Yanggu and HC: Hongcheon) and one southern 
forest (JR: Jirisan National Park) of South Korea (Table 1, 
Fig. 1A). Genomic DNA was extracted from the tissue sam- 
ples by using the DNeasy blood and tissue kit (Qiagen, Val- 
encia, CA, USA), and amplified using the InnoTM Taq 
DNA polymerase kit (Bookyung S-M, Korea). Primers for 
PCR amplification of mitochondrial COI gene was perform- 
ed were Deer F (5'-TTCATTAACCGCTGATTATTTTC 
AAC-3') and Deer R (5’-CACGATATGAGAAATTATACC 
AAACC-3’'). This amplification procedure yielded 680-bp 
fragments of the COI gene; these were aligned using the 
Clustal-X program package (Thompson et al., 1997) and 
manually edited using SeA 1 1.0a (Rambaut, 1996). 
Intraspecific phylogeny was inferred by neighbor-joining 
(NJ) implemented in PAUP* 4.0b10 (Swofford, 2003) and 
Bayesian inference (BI) using MrBayes 3.2 (Ronquist et al., 
2012). The NJ tree was generated using the K2P model, and 
the robustness of the tree was tested with 1,000 pseudorep- 
lications. For the BI analysis, the HK Y model was selected 
as the best-fit model under the Akaike Information Criterion 
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Table 1. List of the samples used in this study 


Intraspecific Phylogeny of the Korean Water Deer 


Common name Scientific name Sampling locality Spree Genbank 
abbreviation accession No. 
Korean water deer Hydropotes inermis argyropus Yanggu, Gangwon Province, South Korea YG1 KP100311 
H. i. argyropus Yanggu, Gangwon Province, South Korea YG3 KP100312 
H. i. argyropus Yanggu, Gangwon Province, South Korea YG4 KP100313 
H. i. argyropus Yanggu, Gangwon Province, South Korea YG5 KP100314 
H. i. argyropus Yanggu, Gangwon Province, South Korea YG6 KP100315 
H. i. argyropus Yanggu, Gangwon Province, South Korea YG7 KP100316 
H. i. argyropus Yanggu, Gangwon Province, South Korea YG8 KP100317 
H. i. argyropus Yanggu, Gangwon Province, South Korea YG9 KP100318 
H. i. argyropus Jiri Mountain Province, South Korea JR16 KP100320 
H. i. argyropus Jiri Mountain Province, South Korea JR17 KP100321 
H. i. argyropus Jiri Mountain Province, South Korea JR18 KP100322 
H. i. argyropus Jiri Mountain Province, South Korea JR19 KP100323 
H. i. argyropus Jiri Mountain Province, South Korea JR20 KP100324 
H. i. argyropus Sokcho, Gangwon Province, South Korea SC15 KP100307 
H. i. argyropus Sokcho, Gangwon Province, South Korea SC21 KP100308 
H. i. argyropus Sokcho, Gangwon Province, South Korea SC22 KP100309 
H. i. argyropus Sokcho, Gangwon Province, South Korea SC24 KP100310 
H. i. argyropus Hongcheon, Gangwon Province, South Korea HC23 KP100319 
Chinese water deer H. i. inermis Jingsu, China JI GQ329030 
H. i. inermis Zhejiang, China ZH1 GQ329027 
H. i. inermis Zhejiang, China ZH2 GQ329028 
Tufted deer? Elaphodus cephalophus Unknown = DQ873526 
Eurasian elk? Alces alces Unknown = EU835716 
@The sequences of Chinese water deer and Eurasian elk were obtained from the GenBank. 
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Fig. 1. A, Collection localities for the Korean water deer. Further information on the samples is provided in Table 1; B, Intraspecific 
phylogeny of Korean and Chinese water deer inferred from the neighbor-joining (NJ) and Bayesian inference (BI) analyses, based 
on mitochondrial cytochrome oxidase I (COI) gene sequences (680 bp). Only the NJ tree is shown in the figure, and the num- 
bers at each node indicate the bootstrap values for the NJ tree (left) and posterior probabilities (shown as percentages) for the BI 
tree (right). Alces alces and Elaphodus cephalophus were used as outgroup species. The Korean water deer formed two separate 
clades, regardless of the collection localities. 
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Clade Subspecies Population 


Nucleotide position Haplotype 
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Fig. 2. Haplotype patterns of Hydropotes inermis. In Korean water deer, two different haplotypes KH1 and KH2 were found; there 
was no sequence difference between KH1 and the Chinese CH1 haplotype. 


Table 2. Pairwise genetic distances between the Korean and Chinese water deer haplotypes 


Chinese water deer 


Korean water deer 


Population Haplotype 
CH1 CH2 KH1 KH2 
Chinese water deer CH1 = > = F 
CH2 0.00466 E = S 
Korean water deer KH1 0.00000 0.00466 = = 
KH2 0.01079 0.01555 0.01191 = 


(AIC) by using Modeltest 3.7 (Posada and Crandall, 1998). 
The BI analysis was performed using the Markov chain 
Monte Carlo technique (MCMC) and one sample every 500 
generations was used. The first 25% from each run was dis- 
carded as burn-in. 


RESULTS AND DISCUSSION 


The NJ and BI analyses generated similar trees with iden- 
tical topologies, with greater support for the NJ tree than 
the BI tree. The analyses confirmed that Korean water deer 
could be separated into two phylogroups (Fig. 1B), one of 
which clustered with the Chinese water deer (JIN, ZH1, and 
ZH2). This result is consistent with that of a previous study 
that analyzed the mitochondrial cytochrome b gene and its 
control region and concluded that some Korean populations 
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might belong to the Chinese water deer (Koh et al., 2009). 
Two haplotypes (KH1 and KH2) were found in Korean wa- 
ter deer (Fig. 2). The nucleotide sequence of the KH1 haplo- 
type did not differ from that of the Chinese CH1 haplotype. 
The K2P distance between the Korean KH1 and KH2 hap 
lotypes was 0.012, and it was considerably higher than the 
distance (0.000-0.005) between the Korean KH1 and two 
Chinese haplotypes (CH1 and CH2) (Table 2). Thus, the 
Korean water deer has at least two genetic lineages, namely, 
H. i. argyropus and another that may belong to the Chinese 
water deer, which was previously thought to be restricted to 
China. The populations of the two Korean genetic lineages 
did not show distinct geographic distributions (Fig. 1A). 

We conclude from our intraspecific phylogeny and haplo- 
type analysis based on the mitochondrial COI gene sequen- 
ce that both the Chinese and Korean water deer are distrib- 
uted across South Korea. Further morphological studies of 
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the Korean water deer will be required to confirm this reas- 
sessment of its taxonomic status. 
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